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ABSTRACT 

Morphologically distinct populations of the Gammarus pulex 
group which should be attributed to at least three species 
applying morphological criteria postulated by Karaman & 
Pinkster (1977) have been reported in the Massif de la 
Sainte Baume (southern France). These populations 
proved to be completely interfertile and thus were con¬ 
sidered conspecific. This group of populations, in which 
three major forms A, B and C were distinguished, has 
been described as Gammarus stupendus Pinkster, 1983. 

Systematic sampling was carried out to determine the 
distribution of the various forms of G. stupendus in the 
Massif de la Sainte Baume and adjacent areas, and the 
relation between morphological and genetic differentia¬ 
tion has been investigated. Genetic differentiation has 
been studied in four polymorphic gene loci: Gpi , Got-1 , 
Mdh-1 , and Mdh-2. 

As the populations used in hybridization experiments 
by Pinkster (1983) proved to be genetically very similar, 
additional crossbreeding experiments have been carried 
out with two genetically more distinct populations. More¬ 
over, experiments were run to test the stability of 
diagnostic characters in each form A, B and C of G. 
stupendus . 

No clear correlation could be found between any of the 
factors morphology, genetic differentiation, and geo¬ 
graphic distribution pattern. In the crossbreeding 
experiments genetically more distinct populations 
appeared to be incompletely interfertile. Thus, the various 
G. stupendus populations, exhibiting different degrees of 
interfertility which may ultimately result in reproductive 
isolation, may indicate different stages in the process of 
speciation. 

RESUME 

Des populations de Gammares du groupe pulex mor- 
phologiquement distinctes qui devraient, en appliquant 
les criteres formulees par Karaman & Pinkster (1977), etre 
attributes au moins a trois especes differentes, ont ete rap- 
portees du Massif de la Sainte Baume (Midi de la France). 
Ces populations, dont trois formes principales A, B et C 
ont ete distinguees et qui ont prouvees d’etre tout a fait 
interfertiles, etaient done considerees conspecifiques et ont 
ete decrites comme Gammarus stupendus Pinkster, 1983. 


Un echantillonage systematique a ete effectue afin de 
determiner la repartition des formes differentes dans le 
Massif de la Sainte Baume et les regions avoissinantes, et 
la relation avec la differenciation genetique. La variation 
genetique de quatre loci polymorphiques a ete etudiee: 
Gpi ; Got-1 , Mdh-1 et Mdh-2. 

Etant donne que les populations utilisees par Pinkster 
(1983) dans des essais d’hybridation etaient tres 
semblables du point de vue genetique, d’autres essais 
d’hybridation ont ete effectues avec des populations 
genetiquement plus differentes. En plus, des essais ont ete 
effectues pour tester la stabilite des caracteres diagnosti- 
ques de chacune des formes A, B et C de G. stupendus. 

Une correlation nette entre ni la differenciation mor- 
phologique et genetique d’une part, ni la differenciation 
morphologique ou genetique et la repartition geographi- 
que d’autre part, n’a pu etre demontree. Dans les essais 
d’hybridation, des populations genetiquement plus dif¬ 
ferentes ne se sont prouvees pas tout a fait interfertiles. 
Done, les populations differentes de G . stupendus , en mon- 
trant plusieurs niveaux d’interfertilite dont pourrait 
ulterieurement resulter une isolation reproductive, 
devraient vraisemblablement etre considerees comme des 
stades differents dans le processus de speciation. 


INTRODUCTION 

Pinkster (1983) investigated morphologically 
different populations of the Gammarus pulex 
group in the Massif de la Sainte Baume and the 
adjacent karst areas in southern France. The 
material he collected could not be identified 
satisfactorily with the aid of the key by 
Karaman & Pinkster (1977). Some populations 
could be identified as G . fossarum Koch, 1835, 
but most of the populations could not be 
attributed to any of the known taxa. Applying 
a set of so-called stable characters (Pinkster, 
1983) at least three species should have been 
distinguished, but even then uncertainty 
remained because of the enormous variability. 
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Although not one population seemed to be 
identical to any other, Pinkster (1983) was able 
to recognize three major forms A, B and C. 
These populations all have certain characters in 
common with G. fossarum, but differ from it 
(and often from each other) in two or more 
characters which should be discriminative 
according to Karaman & Pinkster (1977). Two 
of these forms resemble the French populations 
of G. ibericus Margalef, 1951 (later described as 
G. orinos Pinkster & Scholl, 1984) in more than 
one respect. However, in crossbreeding 
experiments the forms A, B and C proved to be 
interfertile, whereas crosses between any of 
these populations and standard populations of 
“true” G . fossarum and G. orinos failed to give 
any offspring (Pinkster, 1983). Therefore it was 
decided that populations of forms A, B and G 
belonged to a new species, which was described 
as G. stupendus Pinkster, 1983. 

The results of the crosses between the forms 
A, B and C and the morphology of the offspring 
did raise some doubt about the taxonomical 
value of certain morphological characters as 
proposed by Karaman & Pinkster (1977): even 
those which have so far been considered stable 
in some species were too variable in others and 
could not be used for discrimination. 

The primary aim of the present study was to 
determine if morphological differentiation in G. 
stupendus (i.e. the forms A, B and C) coincides 
with genetic differentiation. Second goal was to 
determine if there is any relation between either 
morphological or genetic variation and the 
geographical distribution pattern. 

Genetic variation was investigated by starch- 
gel enzyme electrophoresis. Sampling for elec- 
trophoretical analysis was preceded by 
systematic sampling to map in more detail the 
actual distribution of the morphological forms 
A, B and C of G. stupendus. Thus, sampling car¬ 
ried out in 1980-1981 and earlier (Pinkster, 
1983) was completed and differences in tem¬ 
poral and spatial distribution patterns of these 
forms were checked. 

Laboratory experiments were carried out to 
verify whether the morphological differentia¬ 


tion is constant, i.e. not subject to seasonal 
and/or environmental fluctuations. During 
electrophoretical analysis it was noticed that 
populations used in crossbreeding experiments 
by Pinkster (1983) were genetically rather 
similar to one another. Moreover, “true” G. 
fossarum used in his experiments came from 
northern France and the Alps. In my opinion 
local G. fossarum must also be used in these 
crosses. After the first electrophoretic results it 
was decided to rerun a series of crossbreeding 
experiments with (a) two genetically more 
distinct populations of different morphological 
forms, (b) these populations and a local G. 
fossarum population, and (c) a local G. fossarum 
population and a G. fossarum population from 
northern France. 

STUDY AREA 

As the variability of G. stupendus should be seen 
in the context of a highly structured environ¬ 
ment, the study area will be discussed in more 
detail. The Massif de la Sainte Baume and the 
adjacent limestone areas (fig. 1) include four 
main drainage systems: the Arc, the Argens, 
the Huveaune and the Gapeau. The upper 
reaches of the largest system, the Argens, can 
be subdivided in four smaller tributary systems: 
the Eau Salee, the Cauron, the Caramy and the 
Isolle. Superficial drainages of these systems are 
well separated. However, as some springs in the 
upper reaches are very close to one another (in 
some cases less than 100 m) subsidiary drainage 
into subterranean conduits (a phenomenon 
common in karst areas) through the watershed 
can not be excluded. Thus, an effective 
geographic barrier may be absent. 

The study area is situated in the Mediterra¬ 
nean region and climatical factors as rainfall 
and temperature can be of decisive importance 
to local hydrology. Rainfall is in the order of 
700-1000 mm a year with maxima at the end of 
the winter, in spring and in autumn (Kiener & 
Ollier, 1970; Giudicelli et al., 1980). 
Incidently, rainfall fails in autumn. The area is 
subject to a well-defined “dry season”, ranging 
roughly from June to September. According to 
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Bagnouls & Gaussen (1953) the 4 ‘dry season” is 
defined as the period in which the monthly 
means of rainfall (in mm) is equal to or less than 
twice the temperature (in degrees centrigrade). 
The xerothermic index as defined by these 
authors is 60 for Marseille (for comparison: the 
xerothermic index for Lyon is 0; for Marakech, 
Morocco, it is 178). 

The upper reaches of the rivers Gapeau and 
Argens depend for their water supply directly 
on subterranean groundwater basins and their 
discharge is quite independent on local rainfall. 
The discharge of other tributaries, however, is 
directly related to local rainfall (Kiener & Ollier, 
1970; Giudicelli et al., 1980; pers. obs.). The 
Latay for instance is in the dry season reduced 
to a series of stagnant pools, which often are not 
interconnected anymore, and under extreme 
conditions as prevailing in the summers of 1985 
and 1986 some tributaries desiccate almost 
entirely (e.g. Cauron, Isolle, upper reaches of 
Huveaune and Eau Salee). This is even known 
to happen incidently to the Gapeau system as a 
whole (Kiener & Ollier, 1970). This situation 
can last for several months (up to seven and 
perhaps more; pers. obs. 1985-1986). 

The frequent desiccation of riverbeds in the 
dry season, resulting in the formation of 
isolated river sections and pools subsisting for at 
least several months, may provide a temporal 
genetic barrier. It is obvious however, that 
when complete desiccation occurs, gammarids 
may be absent for quite some time. 

MATERIAL AND METHODS 
Sampling and collection sites 

In November 1984 intensive sampling was carried out in 
the study area to determine the distribution of the various 
morphological forms. In total at 48 localities Gammarus 
were found (table I, fig. 1). In many localities however, no 
gammarids occurred at all. At every locality samples of 
approximately 100 animals have been taken. The animals 
were killed on the spot with 4 % formaldehyde and later on 
transferred to and stored in 70% ethanol. From the 48 
localites, 21 were selected for additional morphological 
and electrophoretic studies. Sampling for this purpose has 
been carried out in November 1985. For comparison two 
populations from the adjacent Massif des Maures were 
sampled as well. 


Animals sampled for electrophoretic experiments were 
anaesthetized with Sandoz MS 222 (1 :1000), identified in 
the field and immediately transferred into a liquid 
nitrogen container (~170°C). In addition to each frozen 
sample, an ethanol control sample was taken for mor¬ 
phological study. As all sampling was done in the winter 
period, in July 1985 additional samples have been taken 
at six localities to check eventual seasonal changes in mor¬ 
phology. 

Material collected by Pinkster (1983) and by J. Berner 
(unpubl.) was reexamined. As far as these localities do not 
correspond with my sampling stations they are indicated 
as stations a-i in table I and fig. 1. 

In November 1984 G. stupendus forms A, B and C of 
representative stations have been collected (stations 1, 32 
and 43) to test the stability of morphological characters 
under uniform laboratory conditions. Part of the animals 
collected were frozen to compare allele frequencies in 
November 1984 and November 1985. 

In November 1985 and July 1986 population samples of 
G. stupendus form B and C, and of G. fossarum were col¬ 
lected at stations 43, 1, and 13, respectively, for additional 
crossbreeding experiments. “True” G. fossarum for these 
experiments were collected in August 1986 in the river 
Slack in Ambleteuse (dept. Pas-de-Calais, northern 
France). 

Laboratory experiments 

After being transported in plastic bags with habitat water 
under oxygen, the animals of each population were kept 
in storage basins of 40 x 50 x 30 cm. Rough stones were 
provided for shelter and air was blown constantly through 
the water. The animals were fed with Stellaria media Vill., 
a common terrestrial weed. 

To test the stability of morphological characters under 
uniform laboratory conditions 25 males and 25 females of 
each population A, B and C of G. stupendus were kept in 
special units. Each unit consisted of two boxes, an outer 
box (25 x 25 x 10 cm) and an inner box (24.5 x 24.5 x 7 
cm). The bottom of the inner box was replaced by a wire 
netting (mesh 750 fim). The adults, unable to pass the 
wire netting, were kept in the inner box. Juveniles sank 
through the meshes and concentrated on the bottom of the 
outer box, so that they could not be eaten by the adults. 
The experiment was done at 15°C and a light/dark period 
of 16/8 hours. Once a month the morphology of 10 
animals of each of the forms A, B and C was examined. 
The same procedure was applied to the FI and F2 when 
raised to maturity. 

The methods used for the crossbreeding experiments 
were identical to those described by Pinkster (1983). The 
B population he used (collected at station 29) was 
genetically very similar to his C population (station 1; 
tables I and II). Therefore, in the present study the station 
29 population has been replaced by the B population from 
station 43. Contrary to Pinkster (1983) an A population 
was not involved. The first series of crossbreeding 




Fig. 1. Study area and collection sites. Stations numbered 1-48 have been sampled for the present study; those numbered a-i have been sampled by J. Berner 

(unpublished). 
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experiments was run in the period November 1985-March 
1986; the second series in the period August 1986-January 
1987. 

Electrophoresis 

Electrophoresis and staining procedures generally fol¬ 
lowed Siegismund et al. (1985). Of 20 enzyme systems 
assayed, four polymorphic enzymes gave good results: 
glucose phosphate isomerase (GPI, E.C. No. 5.3.1.11), 
glutamic oxaloacetic transaminase (GOT-1, E.C. No. 
2.6.1.1), and malate dehydrogenase (MDH-1 and MDH- 
2, E.C. No. 1.1.1.37). The genetic interpretation of the 
variation is inferential. 

The allele frequencies have been used for calculating a 
matrix of Nei’s indices of normalized genetic identity, 1 
(Nei, 1972). Based upon these figures, an UPGMA den¬ 
drogram (Sneath & Sokal, 1973) has been constructed. 

RESULTS 

Distribution of the major morphological forms 
(Table I, figs. 1 and 2) 

The Arc system 

In the Arc system G. stupendus was not found. 
G. pulex gallicus S. Karaman, 1935 dominated, 
and only at three localities (stations 8, 12 and 
13) G. fossarum occurred. As a whole, the Arc 
system is rather polluted (fig. 1) and it is prob¬ 
ably due to this pollution that G. fossarum 
populations regressed and progressively have 
been replaced by populations of G. p. gallicus , a 
species inhabiting generally middle and lower 
reaches, having a higher tolerance to pollution. 

Upper reaches of the Argens system 

In the upper reaches of the Argens system 
(Argens, Eau Salee and Cauron) G. stupendus 
forms A and B, and G. fossarum were found. 
Populations of G. stupendus form A showed only 
to a small extent the diagnostic characters of 
this form as given by Pinkster, 1983 (in par¬ 
ticular the setation of P5-P7) and can therefore 
very easily be confused with G. fossarum (both 
species occur sympatrically at stations 17, 18, 
19, 20, 21 and 22). The morphology of these 
populations was stable and distinct (individuals 
showing transitional characters toward other 
forms of G. stupendus have not been found). 
Thus, at localities where G. stupendus forms A 


and B occur sympatrically (stations 19, 20, 21 
and 22) they could be separated without dif¬ 
ficulties. 

The Caramy subsystem 

In the Caramy G. stupendus form B dominated. 
Although individuals showing transitional char¬ 
acters attributable to both G. stupendus forms A 
and C (as described by Pinkster, 1983) were 
found, no discrete different forms could be dis¬ 
tinguished. 

The Isolle subsystem 

In the Isolle all three forms A, B and C of G. 
stupendus were encountered. Morphological 
characters of G. stupendus form A in the popula¬ 
tion of station 39 were stable, but in both A and 
C populations from station 40 intermediate 
forms showing typical form B characters 
(especially the setation of P5-P7) were fre¬ 
quently found. Material collected by J. Berner 
at the localities c, d and e (see table I and fig. 
1) contained predominantly G. stupendus form A 
identical to the population from station 39. 

The Huveaune system 

The lower reaches of the Huveaune (stations 
34, 35 and 43) are populated by G. stupendus 
form B. The morphology of these populations 
was rather variable. Individuals supposed to 
belong (according to Pinkster, 1983) to G. 
stupendus form A were occasionally found; this 
applies particularly to the populations of sta¬ 
tions 34 and 35. However, as all intermediate 
stages between form B and form A were found, 
this variability can be seen as intra-populational 
variation. The same applies to the G. stupendus 
form A populations in the upper reaches of the 
Huveaune (stations 27, 28, 32 and 33). Thus, 
the G. stupendus form A populations of the upper 
reaches and the form B populations of the lower 
reaches showed an overlap in morphology. 

Older material from the lower Huveaune 
reaches (stations 35 and 43) showed a form A 
morphology and might indicate intermittent 
temporal changes in morphology over longer 
periods as well as shifts in the distribution 
pattern. 
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Table I 


Sampling localities and morphological forms: distribution 
over the various drainage systems (cf. fig. 1 and footnote). 


drainage 

(sub)system 

station 

morphology 
1970- 1980- 

1971* 1981** 

1984 

1985 

Arc 

8 

— 

— 

/ 

— 


16 

— 

— 

p 

— 


9 

— 

— 

p 

— 


14/15 

— 

p 

p 

— 


10/11 

— 

p 

p 

— 


12 

— 

ip 

— 

— 


13 

— 

/ 

/ 

i 

Eau Salee 

24 

— 

B,/ 

/ 

i 


26 

— 

A 

B,/ 

— 


25 

— 

/ 

P 

p./ 


23 

— 

/ 

B 

— 

Argens 

17/18 

— 

A 

A,/ 

A J 


21 

A 

A 

A,/,B 

— 


22 

— 

— 

A,B,/ 

f, B, A 


a 

B 

— 

— 

— 


b 

A 

— 

— 

— 

Cauron 

19/20 

— 

— 

B,A,/ 

A,B,/ 


31 

— 

— 

B ,P 

— 


30 

— 

B 

B ,p 

B 

Caramy 

44 

— 

— 

B 

B 


37 

— 

— 

B 

B 


36 

— 

— 

B 

B 


38 

— 

B 

B, A 

B 

Isolle 

46 

B, A 

B 

O 

B 


39 

— 

— 

A, PJ 

— 


40 

— 

— 

C, A 

C, A 


c 

B 

— 

— 

— 


d 

B, A 

— 

— 

— 


e 

A,B 

— 

— 

— 

Huveaune 

35 

A 

B 

B 

— 


34 

B 

B 

B 

B 


43 

A 

— 

B 

B 


32 

— 

—■ 

A 

A 


33 

— 

— 

A ,p 

— 


27/28 

— 

A 

A 

A 


29 

B 

B 

B 

B 

Gapeau 

3 

— 

— 

C 

C 


2 

— 

— 

C 

C 


1 

C 

C 

C 

G 


4 

C 

— 

C 

— 


45 

A 

— 

C, A 

C,A 


5 

C 

— 

C 

— 


6/7 

A,C 

— 

C, A 

C, A 


41 

B 

— 

A 

O 


42 

A 

— 

A 

— 


f 

A 

— 

— 

— 


g 

A 

— 

— 

— 


h 

A 

A 

— 

— 

— 

Maravenne 

47 

A 

/ 

_ 

B 

Preconil 

48 

— 

B 

— 

B 


Compared to the other sampling sites the 
Fontaine de St. Pilon (station 29) represents a 
particular case. This isolated fountain, situated 
just under the edge of the Sainte Baume moun¬ 
tain crest, apparently does not belong to any 
stream system. However, in periods of rainfall 
water from the fountain basin seeps down 
through the woodlands along the mountain 
slope, most probably toward the upper reaches 
of the Huveaune system. The very small 
population of this locality (“type-locality” of G. 
stupendus form B in Pinkster, 1983) generally 
showed a form B morphology, although 
intergrades with more or less pronounced char¬ 
acters of G . stupendus form A were also found. It 
should be emphasized, however, that some 
diagnostic characters (e.g. setosity of P5-P7) 
provided by Pinkster (1983) revealed to be sex¬ 
ually dimorphic in many populations of G . 
stupendus forms A and B, and their growth 
appears to be allometric (i.e. animals showing a 
form B morphology in the younger life stages 
may end with a form A morphology). For these 
reasons, it was difficult to attribute a certain 
number of individuals to one form or another. 

The Gapeau system 

In the Gapeau system G. stupendus form C 
dominated, although G. stupendus form A was 
also found. The G . stupendus form C popula¬ 
tions showed little variation in morphology. G. 
stupendus form A populations were more 
variable compared to one another. The popula¬ 
tions at stations 6 and 45 exhibited a stable 
morphology, and resembled G. stupendus form 
A populations from the Argens system (stations 
17, 18, 19, 20, 21 and 22). In populations of 
stations 41 and 42 the setosity of P5-P7 was 
more variable, and in this respect these popula¬ 
tions resembled the G. stupendus form A popula¬ 
tions in the upper reaches of the Huveaune (sta¬ 
tions 27, 28, 32 and 33). Older material 
collected at localities f, g, h and i was identical 
to form A populations from stations 6 and 45. 

f=G. fossarum; p = G. pulex gallicus; A = G. stupendus form A; 
B = G. stupendus form B; C = G. stupendus form C; bold 
type = form predominant (figures in order of decreasing 
abundance); O - station desiccated, * = leg. J. Berner; 
leg. Pinkster (1983). 
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Fig. 2. Distribution of the major morphological forms A, B and C of G. stupendus , and G. fossarum. 


The G. stupendus form C populations from the 
Gapeau system were identical to the former 
population from the type-locality of G. 
stupendus , a spring in the Vallat de Fontanieu 
(Pinkster, 1983). This location is nowadays 
canalized and polluted (fig. 1: approx. 10 km 
N.W. of Toulon), and gammarids were not 
found anymore. 

Although G. stupendus form C was largely 
restricted to the Gapeau system whereas the 
tributaries to the Argens were predominantly 
populated by G. stupendus form B populations, 
there does not seem to be a clear correlation 
between the different morphological forms of G. 
stupendus and the drainage systems studied (fig. 


2). Notably the occurrence of G. stupendus form 
A populations across all these systems suggests 
the absence of any correlation. 

Seasonal changes in morphology 

To check seasonal variation in morphology, the 
following stations were sampled: 1 (Gapeau 
resurgence, G. stupendus form C), 17 (Argens 
resurgence, predominantly G . stupendus form 

A) , 27 and 32 (sources of the Huveaune, G. 
stupendus form A), 43 (spring in one of the lower 
reaches of the Huveaune, G. stupendus form B) 
and 29 (Fontaine de St. Pilon, G. stupendus form 

B) . Samples collected in July 1985 at stations 1, 




























Distribution of allelic frequencies at four polymorphic enzyme loci. For legends of morphological forms and station numbers see table I; 
a = fastest moving allele; b-g = slower moving alleles; Sa « silent allele; N = sample size; h - heterozygosity per locus (direct count); H » average 
heterozygosity (direct count) over the four loci studied; (*) = significant departure from Hardy-Weinberg equilibrium (^<0.05). 


8 M. SCHEEPMAKER - DIFFERENTIATION IN GAMMARUS STUPENDUS 




OM 

5 

8 

O'* 

0 



04 

m 

I'D 

8 5 

n 

■CD 

t^i 




CD ^ 



S £ 

-4 

on ^ 

»n 

SP Cft 

i 

*3 

3 

3 

0 

1 

i 

8 3 

1 O 

3 0 

3 1 

3 

3 3 

t 

1 1 

S O 

1 

1 

3 3 1 

3 

S3 3 





£ 

e 





■s 

—r 04 

8 


r^- 

f-j 


a PO 



ro 04 


O * 

m 

'# m 

. 

1 

3 

3 

1 

1 

1 

S3 3 

1 3 

3 3 

1 3 

8 3 

1 3 

3 

1 4 

5 3 

4 

1 

1 3 3 

1 

— 3 



» 

vTr 

CD 

CD 

D- 

3 

t 

S 

3 

1 

1 

on 

— 3 

53 

» 3 

O' 

3 1 

■ 1 

-T 

2 3 

1 « 

4 

1 1 

ro 

— O 

1 

4 

1 -« 1 

4 

m 

CD 

3 

« _> 


■“1 

CD 

3 

SJ 

3 

i 

8 

3 

1 

1 

n 

3 3 

1 1 

04 04 

O O' 

3 3 

3 =» 

3 3 

n 

3 

1 « 

1 

1 1 

O'- 

04 O 

1 

1 

3 8 

3 3 i 

4 

^5 3 

8 

3 




85 





04 

>N 





O' 

m 

■3 


O 



3 


S 

3 


• 

• 

CD 

3 

1 

1 

1 

OO . 

—• O 

1 1 

. - 

1 1 

^ CD 

1 3 

3 

1 i 

£4 3 

4 

1 

3 3 1 

» 

S3 3 

3 



$ 

04 

8 









£ 


O 




? 0 


_ a 


£ 

i 

CD 

CD 

3 

1 

» 

» 

to 0 

1 

• 1 

• « 

to 0 

1 3 

1 

3 1 

to CD 

» 

1 

4 3 3 

4 

*5 3 

3 



3 

# 

a 




04 


> 

_ 

♦ C ' 1 







3^3 


r-, ^ 

£ 

<_> 

» 

CD 

CD 

3 

1 

1 

1 

^ O 

t 1 

1 3 

t 3 

0 - 3 

4 -H 

4 

1 1 

3 CD 

1 

1 

3 3 3 

• 

m 0 

3 



8 

r-4 

O' 

3 




•04 

£ 

3 3 


fl 







3 P5 


m 

<M 

rn 

i 

o 

CD 

3 

1 

1 

» 

* 3 

i 3 

0 3 

1 1 

P3 


4 

1 4 

ro 0 

1 

1 

3 3 i 

4 

$ 3 

3 



o 

£ 


3 



r- 


to 













S aft 

t 

CD 

CD 

i 

3 

1 

1 

04 ^ 

♦ a 

1 3 

1 3 

. . 

^ • 

^ O 


1 

1 1 

t? 0 

1 

1 

1 — 1 

1 

it? 

3 




3 



r- 

0 


■» 

0 


z % 







? a 


.3 


<z 

i 

1 

CD 

i 

1 

3 

3 

S 3 

1 3 

1 3 

■ 1 

3 

1 — 

1 

1 1 

■04 CD 

1 

1 

1 3 3 

4 

3 

3 




CO 

8 





5 

O' 


ni 










_ 

32 

A 

i 

CD 

CD 

3 

1 

1 

i 

It? 3 

1 3 

1 3 

1 • 

ro 3 


1 

1 1 

rn 0 

4 

« 

1 — 1 

4 

01 0 

3 




— 

-4 


0 


* y 


5 m 












3 

m 

i 

CD 

3 

3 

1 

3 

1 

5 3 

0 

1 3 

3 3 

. . 

O^ * 

•T CD 

1 - 

1 

1 1 

8 CD 

1 

4 

4 ^ 4 

4 

5 CD 

3 




£ 

to 

3 

H 


in 

O' 













n 

£ ~ 

i 

« 

3 

3 

3 

3 

1 

S 3 

» 3 

1 3 

t 1 

^ 3 


1 

1 1 

0 

1 

4 

4 — 4 

1 

^ CD 

3 



r~ 

3 

8 

3 



O' 

04 

-n 

04 

rn 

r- 

s 

rn 

CD 

cr 

° 


04 S 


5 

£ 

CD 

^ s 

CD 


» 

3 

3 

3 

3 

1 

1 

♦ 3 

3 » 

1 3 

1 3 

^ 3 

1 3 

3 

4 1 

3 s 3 

1 

3 

1 3 4 

3 

t? 3 

3 





3 

0 

s 

c? 

? 

s 

£ 


3 

8 

5t 

£ 

rn 



S: 3 



5 

3 c 

i 

1 

• 

3 

3 

3 

3 

3? 3 

3 1 

1 3 

1 1 

35 3 

1 3 

3 

* 3 

07 3 

1 

• 

1 3 3 

1 

CD 

3 



3 

85 

8 

3 

3 


PI 

3 *3 

r". ao 

0 

5 8 

£ 

0 


a 

3 



O- ^ 

Sj 

$ _ 

i 

CD 

3 

3 

3 

3 

1 

So 3 

3 3 

3 3 

1 3 

fC 3 

1 3 

3 

1 4 

a 3 

3 

3 

3 3 3 

1 

00 3 

3 



- 

£ 

8 

8 



♦ 

3 

0 0 

S 

3 






04 

CD 

•3 8 


„ m 

s 

3? pet 

i 

CD 

3 

3 

3 

4 

1 

rn 3 

1 3 

3 3 

1 3 

? 3 

» - 

1 

1 4 

S 0 

4 

3 

3 3 4 

4 

3 3 

3 




r-~ 


0 

r^ 

rn 

3 

'3 

rt 

^ 3 


04 

CM 

in 

? 


O' S 

s 


-0 .n 


rn 

_ ♦ 

m 

37 

B 

• 

1 

3 

0 

3 

3 

3 

rn 3 

1 3 

3 3 

« 1 

P5 3 

1 3 

3 

1 1 

St 3 

3 

1 

3 3 1 

• 

o3 3 

3 





lO 




£ £ 





S 

04 

O- 


r»- 
_ m 



3 £ 



8 

I « 

i 

1 

i 

3 

3 

1 

1 

3 

1 t 

1 - 

■ 1 

O O 

1 0 

3 

1 1 

3 

1 

1 

3 3 4 

1 

3 

0 

S in 



CD 

0 

r M 

04 

0-4 


? 

O' 

3 


8 










■04 

i 

1 

3 

3 

3 

3 

1 

35 3 

1 3 

1 3 

1 • 

5 3 

i - 

1 

4 1 

S 0 

4 

4 

4 1 

1 

S 0 

3 




£ 

00 

r— 


rn 


8 

IT 

5 3 


C4 










•n 

o- 

«-h m 

« 

1 

3 

3 

1 

3 

1 

3 

1 3 

3 3 

1 1 

S 3 

1 - 

4 

1 4 

rn 0 

1 

1 

4 — 1 

4 

m a 

3 




o 

3 

co 

3 


0 - 



.? s 

04 










£ 

O' 

— •« 

< 

1 

3 

3 

3 

3 

1 

> • 
in 0 

1 ■ 

4 < 

3 3 

S 3 


1 

1 1 

n cd 

1 

1 

3 3 1 

4 

£ 3 

3 


3- 

~o 

o 






•74 














3 


3 

CD 

i 

1 

1 

t 

1 

00 a 

t 1 

—« 1 

■ 1 

r- O 

1 ~. 

1 

1 4 

!-~ CD 

1 

4 

1 1 

1 

0 

3 

22 

B 



s 

in 

r- 


O' 

0 


£ 5 

cm 


•n 




^ O 





So 

m 

1 

1 

3 

3 

• 

3 

3 

S3 3 

3 3 

1 3 

t 1 

K 3 

1 -H 

' 

1 1 

4 

‘ 

1 — 4 

4 

3 

22 

A 





«n 

n 




in 

in 











3 

1 

1 

i 

1 

3 

3 

1 

♦ — 

• 1 

4 3 

3 1 

m cd 

4 — 

1 

I 4 

n cd 

1 

4 

4 1 

1 

«n 0 

3 


CM 

T“" 


a 







£ 

0 

0 

00 

04 


m 



O ? O 


ro 3 

s 


3 

CD 

3 

1 

i 

1 

» 

8 3 

1 1 

3 1 

3 1 

r® 3 

1 3 

3 

1 < 

3 

4 

< 

3 3 3 

1 

38 3 

3 




o 

£ 

O' 

£ 

0 



- 

£ 2 

m ^ 

S= 

s 


□ 


, 

- 8 


^ S 

1 *— 

17 

A 

1 

1 

3 

3 

3 

3 

3 

” 3 

1 1 

1 3 

3 3 

0 3 

1 3 

3 

4 1 

t3 3 

1 

3 

3 3 * 

1 

<N • 
m cd 

3 

17 

/ 


1 

1 

1 

1 

1 

1 

xi 0 

1 1 

- , 

1 1 

O' 0 

4 w 

4 

4 1 

-0 O 

1 

1 

1 — 1 

4 

-XI 0 

0 


& 

CD 

1 

3 

3 

1 

1 

1 

, 

*8 

i 1 

£ 

3 1 

8 

3 1 

C4 

-a *“! 

— 0 


4 

1 1 

^ <=» 

4 

1 

1 — . 4 

1 


£ 

3 











? 

£ 

0 




8 0 



8 3 


Q a 

3 


— 

1 

i 

• 

1 

1 

1 

3 O 

1 1 

3 1 

3 1 

S 3 

* - 

1 

1 4 

1 

1 

1 3 3 

4 

n 3 

3 











3: 

3 

ao 

3 



in 






3 

13 

t 

- 

1 

> 

• 

t 

• 

i 

1 ^ O 

< 1 

3 1 

3 1 

in * 

■n 0 

4 3 

3 

4 1 

5n 3 

1 

1 

4 — 1 

1 

8 0 

3 


m 

-O 

u 

■0 

01 


O* 

m. 4 = 

<* jo 

id a 

ai 

* 

f» a 

u 

T3 • 

* ^ 

<• 

•a 

nut> 

9* 

X 4= 

X 

S A 









T-1 




r^i 




<M 






Z £ 

Z o 

w e 

<§• 







<s 




i 




i 
















BIJDRAGEN TOT DE DIERKUNDE, 57 ( 1 ) - 1987 


9 


17, 29 and 43 showed no differences compared 
to those collected in November 1984 and 1985; 
samples from stations 27 and 32 showed some 
minor differences in the setation of P5-P7. 

Stability of morphological characters under uniform 
laboratory conditions 

The animals for this experiment were collected 
at stations 1 (G. stupendus form C), 32 (G. 
stupendus form A) and 43 (G. stupendus form B). 
The experiments were run from November 
1984 to August 1985. The morphology of all 
three forms proved to be stable. Neither the 
animals collected, nor their progeny in FI and 
F2 did show any deviation with regard to the 
initial morphology. However, the setosity of 
P5-P7 of males in form A and B populations 
diminished as they grew older. This phe¬ 
nomenon is much less expressed in females, 
especially those of form B. 

Genetic variation 

Variation in allelic distribution pattern 

All enzymes, MDH-2 excepted, moved 
anodally. Major differences in allele distribu¬ 
tion occurred at the Gpi and Got-1 loci (table II). 
The allele distribution in G. fossarum at these 
loci differed almost completely from all the 
other populations studied; thus, these loci can 
be considered genetic markers for G. fossarum 
and G. stupendus. 

Morphological forms of G. stupendus occur¬ 
ring sympatrically showed different allele 
distribution patterns (table II: stations 6, 19, 
22, 40 and 45). Moreover, it can be seen that 
intra-areal populations do not necessarily exhibit 
a more similar allele distribution pattern than 
inter-areal populations (as reported by Gooch & 
Hetrick, 1979), and only part of the discon¬ 
tinuities in the distribution of allele frequencies 
coincides with differences in morphology. 

Variation among morphological forms 

Values of I (Nei’s index of normalized genetic 
identity; Nei, 1972) range from 1, indicating 
identical allele frequencies, to 0.12, meaning 


that the complements of alleles only slightly 
overlap (table III). As in these indices identical 
monomorphic loci are not included in the 
averages (which increases the discriminatory 
power), these figures are not directly com¬ 
parable with the ranges of indices characteristic 
of specific or subspecific differentiation (Ayala, 
1975; Thorpe, 1982). 

If morphological differentiation in the 
various forms of G. stupendus (viz. forms A, B 
and C) is paralleled by genetic differentiation, 
intra-form values of I should be higher than 
inter-form values. Following Gooch & Hetrick 
(1979) mean / and standard deviation of the 
genetic indices of each morphological form A, B 
and C of G. stupendus , and G. fossarum were 
calculated. These figures and the range lowest- 
highest value, Imin'Imax ( to visualize the 
heterogeneity of 7), are given in table IV. G. 
fossarum , showing the highest mean intra-form 
value of 7, seems to represent a well-defined, 
genetically coherent group of populations in the 
area under study. 

Differences in mean intra-form values of 7 in 
the three forms A, B and C of G. stupendus (table 
IV) are due on the one hand to the higher level 
of genetic differentiation among the popula¬ 
tions showing a form A morphology (resulting 
in low mean intra-form I values), and on the 
other hand to the effect of high I values of 
neighbouring form B and C populations. 
Although minimum-maximum ranges of I in 
form B and C populations are still large (table 
IV), mean intra-form values of I can still be 
high as a consequence of a greater number of 
single high I values of intra-areal comparisons. 
However, as mentioned above, genetically dif¬ 
ferent populations may occur sympatrically 
(tables I and II), and intra-area I values do not 
necessarily have to be higher. 

Compared to G. fossarum , the mean intra¬ 
form value of I in form A of G. stupendus was 
low, with a wide range of I m i n -Imax (table IV). 
Moreover, the mean inter-form A-C value of I 
is higher than the intra-form A value (table IV). 
Using the Mann-Whitney U-test, intra-7 values 
of form A do not differ significantly at the 5% 
level from inter-A-B values (P = 0.71) and inter- 
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Table III 

Matrix of Nei’s indices of normalized genetic identity (/), based upon allele frequencies at four polymorphic enzyme 
loci. Stations in the same order as in table II, legends of morphological forms as in table I. 


station 13 24 25 17 17 22 22 22 19 19 19 30 44 37 38 46 40 40 3* 43 32 27 29 3 1 45 45 6 6 47 48 

■orph. /ABBBBBBACBBAAB C C A C C A B B 


13/ 

24/ 0.76- 

25 / 0.76 1 - 

17/ 0.76 1 1 - 

17 A 0.29 0.54 0.56 0.53 - 

22 f 0.82 0.97 0.97 0.97 0.54 - 

22 A 0.30 0.59 0.6 0.53 0.96 0.53 - 

22 B 0.28 0.55 0.56 0.55 0.6 0.55 0.68 - 

19/ 0.76 1 1 1 0.54 0.98 0.59 0.56 - 

19 A 0.27 0.53 0.54 0.52 0.99 0.51 0.88 0.57 0.52 - 

19 B 0.3 0.56 0.57 0.56 0.59 0.56 0.65 1 0.57 0.56 - 

30 B 0.28 0.54 0.55 0.54 0.64 0.54 0.68 0.99 0.55 0.62 0.99 - 

44 B 0.49 0.36 0.37 0.36 0.57 0.43 0.66 0.56 0.37 0.5 0.53 0.53 - 

37 B 0.43 0.44 0.45 0.44 0.51 0.49 0.64 0.89 0.45 0.47 0.88 0.67 0.59 - 

38 B 0.12 0.4 0.42 0.4 0.45 0.4 0.48 0.7 0.41 0.43 0.72 0.72 0.26 0.7 - 

46 B 0.3 0.58 0.59 0.58 0.58 0.58 0.64 0.79 0.59 0.55 0.8 078 0.44 0.75 0.89 - 

40 A 0.52 0.28 0.29 0.20 0.54 0.37 0.6 0.41 0.29 0.48 0.37 0.4l 0.96 0.52 0.17 0.34 - 

40 C 0.29 0.55 0.56 0.54 0.62 0.55 0.75 0.66 0.55 0.57 0.6* 0.61 0.65 0.63 0.66 0.87 0.55 - 

34 B 0.28 0.55 0.56 0.55 0.6 0.55 0.69 1 0.56 0.58 1 0.99 0.54 0.9 0.73 0.01 0.4 0.66 - 

43 B 0.33 0.62 0.63 0.61 0.65 0.62 0.76 0.97 0.62 0.66 0.96 0.95 0.62 0.88 0.76 0.89 0.47 G.Bl 0.97 - 

32 A 0.27 0.52 0.54 0.52 0.59 0.53 0.74 0.65 0.53 0.54 0.63 0.6 0.67 0.63 0.66 0.B6 0.56 0.99 0.65 0.81 - 

27 A 0.29 0.56 0.57 0.56 0.62 0.56 0.77 0.76 0.56 0.58 0.7* 0.73 0.6 0.74 0.78 0.94 0.5 0.96 0.77 0.89 0.96 - 

29 B 0.28 0.54 0.55 0.53 0.6 0.54 0.72 0.7 0.54 0.56 0.6* 0.67 0.6 0.68 0.76 0.93 0.5 0.97 0.71 0.85 0.98 0.99 - 

3 C 0.12 0.4 0.41 0.4 0.42 0.39 0.46 0.67 0.4 0.41 0.68 0.68 0.25 0.68 0.99 0.9 0.16 0.69 0.7 0.74 0.69 0.8 0.79- 

1 C 0.27 0.52 0.54 0.52 0.59 0.52 0.73 0.62 0.53 0.54 0.6 0.57 0.65 0.6 0.66 0.85 0.55 0.99 0.62 0.79 1 0.96 D.98 0.69 - 

45 A 0.29 0.51 0.52 0.51 0.55 0.57 0.59 0.84 0.53 0.53 0.84 0.86 0,4 0.77 0.77 0.82 0.32 0.64 0.86 0.85 0.64 Q.7d 0.72 0.74 0.61 - 

45 C 0.27 0.51 0.52 0.51 0.58 0.51 0.73 0.62 0.51 0.53 0.6 0.58 0.69 0.62 0.64 0.83 0.58 0.98 0.63 0.79 1 0.95 0.97 0.67 1 0.59 - 

6 C 0.27 0.52 D.54 0.52 0.6 0.52 0.7* 0.61 0.53 0.55 0.59 0.57 0.65 0.61 0.68 0.85 0.55 0.98 0.61 0.78 C.99 0.95 0.97 0.71 0.99 0.6 0.99 - 

6 A 0.33 0.59 0.6 0.58 0.58 0.6 0.62 0.79 0.59 0.56 0.81 0.81 0.39 0.75 0.86 0.9B 0.31 0.81 0.82 0.8S U.79 0.91 0.89 0.88 0.79 0.86 0.76 0.78 - 

47 B 0.38 0.49 0.51 0.49 0.51 0.52 0.57 0.72 0.5 0.48 0.74 0.72 0.5 0.78 0.83 0.96 0.44 0.83 0.75 0.83 0.82 0.91 0.91 0.85 0.82 0.78 0.81 0.81 0.95 - 

48 B 0.25 0.33 0.34 0.33 0.36 0.36 0.39 0.61 0.34 0.34 0.63 0.63 0.32 0.65 0.7 0.70 0.28 0.6 0.64 0.67 0.6 0.71 0.7 0.72 0.59 0.72 0.57 0.59 0.81 0.81 


Table IV 


Mean intra- and inter-morphological form values of the genetic identity (!) and standard deviation, based upon four 
polymorphic enzyme loci. Intra-form values in bold face; range l m in^max * n italics. 



G. stupendus form A 

G. stupendus form B 

G. stupendus form C 

G. fossarum 

G. stupendus form A 

G. stupendus form B 

G. stupendus form C 

0.65 ±0.18 

(0.31-0.99) 

0.65 ±0.18 
(0.17-0.99) 

0.69 ±0.19 
(0.16-1) 

0.76 ±0.16 

(0.28-1) 

0.7 ± 0.14 
(0.25-0.99) 

0.87 ± 0.16 

(0.67-1) 


G. fossarum 

0.49 ±0.17 
(0.27-0.6) 

0.47 ±0.11 
(0.12-0.63) 

0.46 ±0.11 
(0.12-0.56) 

0.92 ±0.1 

(0.76-1) 


A-C (/ > =0.40) values of /. Thus, the mor¬ 
phological form A is not paralleled by a specific 
genetic differentiation. 

The mean intra-form value of / in G. 
stupendus form B is higher than in form A, but 
still lower than in G. fossarum and single I values 
are very heterogeneous. Intra-form values of / 


in form B are not significantly different from 
the inter-B-C values at the 1 % level 
(P= 0.0174), but differ significantly from the 
inter-A-B values (P= 0.0006). This means that 
part of the form B populations shows a uniform 
genetic differentiation. 

The mean intra-form value of I in G . 
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Table V 


Mean intra- and inter-zone values of the genetic identity (7) and standard deviation of four zones of high genetic iden¬ 
tity, based upon four polymorphic enzyme loci. Intra-zone values in bold face; range I m i n -l m ax i n italics. 


zone 

Huveaune-Gapeau 

Huveaune-Cauron- 

G. stupendus form A 

G. fossarum 



Argens 

Argens 

Argens 

Huveaune-Gapeau 

0.98 ±0.02 

(0.95-1) 




Huveaune-Cauron- 

0.68 ±0.09 

0.98 ±0.02 



Argens 

(0.56-0.89) 

(0.95-1) 



G. stupendus form A 

0.63 ±0.08 

0.63 ±0.05 

0.99 ±0.01 


Argens 

(0.53-0.77) 

(0.56-0.76) 

(0.99-1) 


G. fossarum Argens 

0.54 ± 0.18 

0.57 ±0.03 

0.55 ±0.03 

0.99 ± 0.01 


(0.5-0.57) 

(0.54-0.65) 

(0.51-0.6) 

(0.99-1) 


stupendus form G is relatively high. Intra-form 
values differ significantly from inter-A-C 
(P= 0.009) and inter-B-C values (/*= 0.001) at 
the 1 % level, which means that most of the 
form G populations show a uniform genetic dif¬ 
ferentiation. It should be kept in mind, how¬ 
ever, that different morphological forms may 
also be subject to similar genetic differentiation 
(see table II). 

Mean inter-form values of the three forms A, 
B and C of G. stupendus , and G. fossarum (viz. I 
inter-A- G. f , B-G. / and C-G. /.) do not differ 
significantly at the 5% level ( J P= 0.29, 0.12 and 
0.013, respectively), which means that on the 
average all three forms A, B and G of G. 
stupendus differ genetically to about the same 
extent from G. fossarum. 

Genetic variation among areas 

To examine in more detail the distribution of 
populations having a high degree of genetic 
identity, these populations were pooled and 
grouped into “zones of high identity’’ (fig. 3), 
and mean I values of these zones were 
calculated (table V). Since genetically different 
forms of G. stupendus , and G. fossarum occurred 
sympatrically, these zones may overlap partly 
or even entirely (fig. 3). The various zones and 
the morphological forms they include are listed 
below. 

Huveaune-Gapeau: stations 27 and 32 ( G. stupendus form 
A); station 29 (form B); stations 1, 6, 
40 and 45 (form C). 


Huveaune-Cauron- 

Argens: stations 19, 22, 30, 34 and 43 (all G. 
stupendus form B). 

G. fossarum Argens: stations 19, 22, 24 and 25. 

G. stupendus form A Argens: stations 17, 19 and 22. 

It is striking that populations of different mor¬ 
phological forms occurring in the same stream 
system never belonged to the same zone of high 
identity (applies to all zones). To put it other¬ 
wise, populations of different morphological 
forms belonging to the same zone of high iden¬ 
tity never occurred in the same stream system. 
It should be noticed as well that the extension 
of these zones may range across several stream 
divides. High genetic similarity across stream 
divides was also observed in populations of the 
neighbouring upper reaches of the Latay, 
Caramy (stations 3 and 38; / 3 , 3 s = 0.99) and 
I soli e (station 46; ^45 = 0.9, 738,46 = 0.89). 
The same was observed in the Gapeau and 
Isolle (stations 45 and 40; 745,40 = 0.98). 

When genetically different populations occur 
sympatrically, it is evident that discontinuities 
in allele distribution and thus lower 7 values will 
be found within the same stream system. How¬ 
ever, dissimilarities in allele distribution occur¬ 
red also in the same stream system among 
genetically less differentiated populations not 
susceptible to be attributed to different 
genetical or morphological forms (viz. stations 
37 and 38: 737 3 $ = 0.69; stations 1 and 3: 
I \,3 = 0.69; station 40 population C and station 
46: 740,46 = 0.87). Thus, 7 values among 
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Fig. 3. Distribution of four zones of high genetic identity. The zones are composed of 
populations with high mutual genetic identity (I) levels, based upon four polymorphic 
enzyme loci. G. / Arg. = G. fossarum Argens; G. stupendus f. A. Arg. = G. stupendus form A, 
Argens; Arg.-Huv. - Argens-Cauron-Huveaune; Huv.-Gap. = Huveaune-Gapeau. 


populations across a stream divide may be 
higher than among populations in one and the 
same tributary. 

The UPGMA dendrogram 

The results discussed above are summarized in 
an UPGMA dendrogram (cf. Sneath & Sokal, 
1973) (fig. 4). Two main clusters (G. stupendus 
and G. fossarum ) can be distinguished. The G. 
stupendus cluster can be subdivided in smaller 
clusters. The smallest cluster consists of popula¬ 
tions of stations 40 (form A) and 44. These 
populations are somewhat more closely related 


to G. stupendus than to G . fossarum and may have 
more affinities to populations that have not 
been investigated electrophoretically (e.g. sta¬ 
tions 39, 41 and 42, localities a-i; table I, fig. 1). 
The same applies to the population of station 
48, from the Massif des Maures. Other clusters 
cover the zones G . stupendus form A Argens, 
Huveaune-Cauron-Argens and Gapeau- 
Huveaune (table V, fig. 4). 

A number of populations, that previously 
could not be attributed to a zone of high iden¬ 
tity because of relatively low I values with 
respect to all other populations studied, are 
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A: G. stupendus f.A 
B: G. stupendus f.B 
C: G. stupendus f.C 
G f : G. fossarum 


d 



c 


22 

B 

34 

B 

19 

B 

30 

B 

43 

B 

37 

B 

45 

A 

38 

B 

3 

C 

46 

B 

6 

A 

47 

B 

40 

C 

32 

A 

45 

C 

1 

c 

6 

c 

27 

A 

29 

B 

48 

B 

17 

A 

19 

A 

22 

A 

44 

B 

40 

A 

13 

G. 

24 

G. 

25 

G. 

17 

G. 

19 

G. 

22 

G. 


0.4 


0.5 


0.6 


0.7 

I 


0.8 


09 


10 


Fig. 4. UPGMA dendrogram based upon four polymorphic enzyme loci. I = Nei’s index of 
normalized genetic identity; for legends of station numbers see table I and fig. 1. 


included in the dendrogram. However, as the 
cophenetic coefficient is relatively low (0.864), 
their position is arbitrary. For this reason, and 
because of sampling errors due to small sample 
sizes, the dendrogram becomes less reliable 
close to its extremities. However, from this den¬ 
drogram it becomes clear that there is no cor¬ 
relation between morphological and genetic 
variation, and that there is no distinct relation 
between either morphological or genetic varia¬ 
tion and the stream systems studied. 

Comparison of allele frequencies in 
1984 and 1985 

Allele frequencies at the Gpi and Got-1 locus in 
populations of stations 1, 32 and 43 are 
presented in table VI. P values indicate that the 
allele frequencies at these loci in November 
1984 and November 1985 did not differ 
significantly in any of the three populations 
studied. 


Cross-breeding experiments 

Local G. fossarum and G. fossarum from northern 
France (// and f n in tables VI and VII) proved 
to be completely interfertile in reciprocal 
crosses: the F2 (FI x FI) could be obtained. In 
crosses between local G. fossarum and G . 
stupendus form C only one female gave offspring 
(cross C x//, table VIII). From crosses between 
local G. fossarum and G. stupendus form B several 
females with offspring could be obtained 
(fix B); in the reciprocal cross only one female 
gave offspring. FI from these crosses (except for 
C x//) have been raised to maturity, but failed 
to give offspring. Experiments to hybridize G. 
stupendus form B and C were more successful. 
The most striking result from these experiments 
is that although several dozens of FI juveniles 
of both crosses B x C and C x B have been 
raised to maturity, no F2 (FI x FI) could be 
obtained. 
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Table VI 


Allele frequencies at the Gpi and Got-1 locus in populations of stations 32, 43 and 1 in November 1984 and November 
1985; b = fastest moving allele; N « sample size; P » probability in Chi-square contingency test (* = probability in 
Fisher’s Exact test). 





Gpi 





Got-1 




b 

c 

d 

f 

N 

P 

b 

c 

d 

f 

N 

P 

station 32 (1984) 

0.02 

0.88 

0.1 — 

48 1 

| 0.043 

0.04 

— 

0.96 

— 

48 1 

| 0.5* 

(1985) 

0.1 

0.82 

0.08 — 

33 J 

0.06 

— 

0.94 

— 

33 j 

station 43 (1984) 

— 

0.18 

0.61 

0.21 

31 1 

0.53 

0.62 

0.03 

0.35 

— 

16 1 

| 0.6 

(1985) 

0.01 

0.31 

0.52 

0.16 

44 J 

0.62 

0.08 

0.3 

— 

42 J 

station 1 (1984) 

0.04 

0.84 

0.12 — 

34 1 

f 0.16 

— 

— 

0.92 

0.08 

69 1 

| 0.92 

(1985) 

0.03 

0.94 

0.03 — 

74 j 

— 

— 

0.9 

0.1 

74 J 


DISCUSSION AND CONCLUSIONS 

The situation described in this study resembles 
in more than one respect the findings of Gooch 
& Hetrick (1979), who studied the genetic struc¬ 
ture of G. minus in a North American karst 
area. These authors state that genetic similarity 
may be largely a neighbourhood effect quite 
independent of ecophenotypic differentiation 
(in this study: population 1 of form A is likely 
to have a higher / with population 2 of form C 
in the same area than with population 3 of form 
A in another area); they conclude that the forms 
are independently evolved phenotypes at most 
localities rather than genetically knit distinct 
“varieties”, subspecies or species. In their 
opinion it is shown as well that convergent mor¬ 
phologies can evolve without an accompanying 
restructuring of the entire genome toward a 
greater similarity. However, based upon three 
loci, their statement might be somewhat 
premature. Nevertheless, the results of the 
present study (based upon four loci) appear to 
indicate the occurrence of morphological dif¬ 
ferentiation without restructuring the entire 
genome (e.g. populations of G. stupendus forms 
A, B and C, placed in the same cluster in fig. 
4 and in the same zone of high identity in table 
V). These results indicate as well that con¬ 
vergent morphologies may evolve in¬ 
dependently, or that genetic restructuring to 
some degree may occur without accompanying 
changes in morphology (e.g. G. stupendus form 


B populations in Cauron, Caramy, Isolle and 
Preconil; the various form A populations in 
Argens, Cauron, Huveaune, Isolle and 
Gapeau). 

Otherwise, the results of the present study 
are more complex than the situation reported 
by Gooch & Hetrick (1979) and there is one 
essential difference: in the area under study 
genetically and morphologically different forms 
occurred sympatrically. Different mor¬ 
phological forms occurring (whether or not 
sympatrically) in the same stream system were 
all genetically different (e.g. G. stupendus forms 
A and B in the Argens) and sometimes popula¬ 
tions within the same stream system attributed 
to the same morphological form differed also 
genetically (e.g. G. stupendus form B in the 
Caramy). In this context for instance, it is likely 
that G. stupendus form A from the Argens- 
Cauron represents such a genetically-knit 
distinct “variety” (as mentioned above) in 
which morphological changes have been accom¬ 
panied by genetic differentiation, or vice versa. 
Occurring in sympatry with genetically dif¬ 
ferent populations of G . stupendus form B (sta¬ 
tions 19 and 22), these populations should prob¬ 
ably be considered as a separate species. 
However, there is no certainty that the mor¬ 
phology is a direct expression of the changing 
genome. Therefore, based upon the four loci 
studied, it can be concluded that there is no 
direct relationship between morphological and 
genetical differentiation. 
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Table VII 

First series of crossbreeding experiments: November 
1985-March 1986; // = local G. jossarum population from 
station 13; / n = “true” G. Jossarum population from 
Ambleteuse (northern France); B * G. stupendus form B 


(station 43); C = 

G. stupendus form C (station 1). 

cross 

ax 9 

N 

aa 

N 

99 

N 

ovig. 

99 

N 99 

with 

offspring 

F2 

obtained 

fl*fl 

30// 

36/, 

16 

8 

+ 

BxB 

37 B 

40 B 

13 

3 

+ 

C x C 

34 C 

34 C 

28 

14 

+ 

ft*fn 

42 ft 

42/„ 

35 

12 

+ 

fnx.fi 

48 /„ 

50/, 

28 

4 

+ 

fl* B 

35 ft 

22 B 

10 

2 

— 

B x// 

62 B 

59/, 

9 

1 

— 

//xC 

50// 

50 C 

25 

— 

— 

Cx/ ; 

52 C 

52/, 

7 

— 

— 

BxC 

77 B 

77 C 

13 

2 

— 

CxB 

78 B 

65 B 

10 

1 

— 


Table VIII 

Second series of crossbreeding experiments. For legends 
see table VII. 

cross 

a x 9 

N 

aa 

N 

99 

N 

ovig. 

99 

N 99 

with 

offspring 

F2 

obtained 

//X/, 

30/, 

45/, 

14 

6 

+ 

fn x fn 

32 /„ 

32/ B 

15 

10 

+ 

BxB 

98 B 

128 B 

117 

20 

+ 

CxC 

78 C 

106 C 

78 

7 

+ 

II x fn 

50/, 

80/„ 

17 

7 

+ 

fn x fl 

100/„ 

100/, 

34 

9 

+ 

Z/xB 

100// 

108 B 

37 

3 

— 

B x/, 

94 B 

126// 

38 

— 

— 

/,xC 

100// 

117 C 

28 

— 

— 

C x/, 

54 C 

73// 

12 

1 

— 

BxC 

50 B 

68 C 

28 

14 

— 

CxB 

100 C 

100 B 

43 

15 

— 


A clear correlation between morphological 
variation and geographical distribution pattern 
could not be found (fig. 2). The genetical varia¬ 
tion pattern is more coherent (fig. 3) and cor¬ 
responds fairly well to the main clusters of the 
UPGMA dendrogram (fig. 4). The validity of 
the “zones” of high identity and the main 
clusters of the UPGMA dendrogram are con¬ 


firmed by genetic differentiation at Ldh and Pep 
loci (Scheepmaker, in prep.). However, discon¬ 
tinuities in the allelic distribution do not coin¬ 
cide with apparent geographic barriers. 

A number of populations had low /-levels 
with all other populations studied and therefore 
did not figure in any of the zones of high iden¬ 
tity: the G. stupendus form A population from 
the Isolle (station 40) and the form B popula¬ 
tions from Caramy (station 44) and Preconil 
(station 48, Massif des Maures). Possibly these 
populations are more closely related to the 
populations which have not been investigated 
electrophoretically (e.g. populations at stations 
39, 41, 42 and at localities a-i; cf. table I and 
fig. i). 

Unfortunately G. stupendus form A from the 
Argens and the Cauron was not involved in the 
hybridization experiments. One of the main 
reasons for this was that this form always occurs 
sympatrically with G. Jossarum , whereas it is 
very difficult to separate the two forms in the 
field. The taxonomic status of these populations 
will be treated in a future publication (Scheep¬ 
maker, in prep.). 

Hybridization experiments show a close rela¬ 
tionship between the forms considered in this 
study: all crosses gave offspring in at least one 
of the reciprocal combinations. However, 
crosses between G. Jossarum and G. stupendus 
populations (stations 13, 43 and 1; 

^13,43 = 0.33, ^ 1,13 = 0.27) were not very suc¬ 
cessful, and F2 could never be obtained 
(whether or not due to an artefact as a conse¬ 
quence of the limited number of FI individuals 
that could be raised to maturity) and reproduc¬ 
tive isolation between these species has most 
probably been achieved. 

Results of the hybridization experiments 
improve stepwise when inter-populational 
crosses of G. stupendus are considered. 
Reciprocal crosses between form B and C (sta¬ 
tions 1 and 43; 7^43 = 0.79) gave a numerous 
offspring, which could be raised to maturity. 
However, an F2 could not be obtained. Thus, 
unless these results are due to laboratory 
artefacts, reproductive isolation has probably 
been achieved. Finally, in hybridization 
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experiments carried out by Pinkster (1983), the 
offspring (FI) from all crosses between G. 
stupendus form A, B and C (collected at stations 
1, 27 and 29; 7 i j27 = 0.96, 7i, 29 = 0.98, 

^ 27,29 = 0.99) were completely interfertile. The 
results of these experiments suggest several 
genome differences leading toward reproduc¬ 
tive isolation among populations of G. stupendus. 

According to Gooch & Hetrick (1979), high 7 
values may be due to high gene flow, concor¬ 
dance of selective forces or, unless dealing with 
numerous loci, simple coincidence, whereas a 
low 7, even based on a single locus, means 
incontrovertibly that the populations are 
isolated or experience different modes of selec¬ 
tion. Unless coincidence, inter-areal popula¬ 
tions (e.g. forms A and C, stations 27 and 1) 
separated by an obvious geographic barrier 
(over a long distance, across a stream divide) 
may have high 7 values and be completely 
interfertile, whereas other populations (e.g. 
those at stations 1 and 43), apparently 
separated by one and the same barrier do not 
seem to be completely interfertile anymore, 
thus advancing toward reproductive isolation. 

On the one hand, the high /-level of stations 
3 and 38 (7=0.99) indicates that, unless coin¬ 
cidence, gene flow across stream divides is quite 
well possible. Although the upper reaches are 
very near to each other, they desiccate almost 
completely in the dry season. It is not known 
whether these reaches are recolonized during 
the “rainy season” or not; otherwise gene flow 
must occur by subterranean contact. Evidence 
is against a fast recolonization of dried-up water 
reaches: many localities where no gammarids 
were found in 1980, 1982 and 1984 (the “zero” 
samples in fig. 1) appeared to desiccate 
periodically in 1985. 

On the other hand, populations in the same 
tributary showed discontinuities in the allelic 
distribution pattern. 7 values may be lower 
while no obvious genetic barrier is present (e.g. 
populations at stations 37 and 38, / 37 . 3 s = 0.69; 
at stations 1 and 3, 71,3 = 0.69). These discon¬ 
tinuities seem to be correlated with differences 
in average heterozygosity over the four loci 
studied. Average heterozygosity (H in table II) 


among populations in the same zone of high 
genetical identity is generally rather constant. 
Partial reproductive isolation might explain the 
maintenance of different allele frequencies in 
the same system where no obvious genetic bar¬ 
rier is present. 

It can be seen from table I that changes (how¬ 
ever slight) occur in the temporal distribution 
pattern of the morphological forms (e.g. sta¬ 
tions 26, 23, 21, 46, 35, 45 and 47). These 
changes may represent simple shifts in relative 
abundance of the various forms due to chang¬ 
ing environmental conditions favouring one 
form at the expense of another, or the coin¬ 
cidental replacement of one form by another 
when a tributary is recolonized, or both at the 
same time. Thus, the distribution of the various 
forms can be seen as a dynamic process, gov¬ 
erned by intermittent environmental factors. 
For instance, many populations are repeatedly 
separated during longer periods of drought, 
breaking down migration routes and gene flow. 
As has been demonstrated by Goedmakers 
(1981) and Goedmakers & Pinkster (1981), 
gammarids are actively and sometimes 
massively migrating organisms. In northern 
France, migration maxima occur in summer; in 
winter the migration almost stops. If, in the 
area under study, migration maxima occur in 
late spring, summer or autumn, they will often 
coincide with periods of drought. Life cycles in 
lower reaches are possibly adapted to local 
environmental conditions, but in the upper 
reaches or in a source or cave environment, 
where temperature is more constant, 
temperature fluctuations may approximate 
conditions found in similar habitats in central 
and northern France. 

Putting together all the data presented it can 
be concluded that the various forms of G. 
stupendus , showing different degrees of interfer¬ 
tility which may ultimately result in reproduc¬ 
tive isolation, must probably be seen as dif¬ 
ferent stages in the process of speciation. 

In combination with isolation, genetic drift is 
often referred to as a mechanism to provoke 
discontinuities in the allele distribution pattern 
(Bulnheim, 1985; Bulnheim & Scholl, 1981a & 



BIJDRAGEN TOT DE DIERKUNDE, 57 ( 1 ) - 1987 


17 


b, 1986; Siegismund et al., 1985). However, 
unless in the case of extreme bottlenecking, the 
effect of genetic drift is not sufficient to explain 
the differences observed. According to Ayala 
(1982), gene frequency changes will be gov¬ 
erned primarily by random drift if 4 Nx< 1 
(N '= population size; x = rate of other deter¬ 
ministic processes). For example, if we assume 
a migration rate (m) of 2 % (or two individuals 
for every 100) and there is no mutation or selec¬ 
tion, gene frequencies will change toward the 
frequencies in the population from which the 
migrants came, even in a small population with 
only 25 individuals, because in such a case 
4Nm = 4x25x2% = 2>1. Computer simula¬ 
tion experiments by Bodmer & Cavalli Sforza 
(as discussed by Goodenough, 1984), with an 
hypothetical population of 25 individuals and 
two alleles with frequency 0.5, demonstrated 
that in the absence of other deterministic pro¬ 
cesses it may require 42 generations to get one 
of the alleles fixed. Similar experiments have 
been reported by Spiess (1977): 50 populations 
of 20 individuals and two alleles with frequency 
0.5 came all to fixation by no later than the 
94th generation and no earlier than the 6th. In 
the case of G. stupendus this represents respec¬ 
tively 17-24, and 1.5-55 years (Scheepmaker, 
unpublished data). The Fontaine de St. Pilon 
(station 29) for instance, inhabited by an 
isolated population of several dozens to several 
hundreds of individuals (depending on the 
period of the year) has still an average 
heterozygosity of 0.1 (table II), which is not 
exceptional compared to other populations. 

However, the inter-area populations having 
high /-levels and being interfertile are all large 
and inhabit a stable environment (i.e. sources 
and upper reaches with a rather constant 
discharge). Thus, they are not subjected to 
important bottlenecking and less susceptible to 
the effects of genetic drift, whereas selective 
modes may be quite similar. 

Otherwise, local conditions as substrate and 
temperature regime may be very variable (e.g. 
muddy and polluted sections in the lower 
reaches or muddy stagnant pools in desiccating 
riverbeds, in contrast to the upper reaches and 


springs with a constant discharge, stony 
substrate and constant temperature). 
Temperature in or in the neighbourhood of 
springs and caves will be rather constant (yearly 
fluctuations in the order of 1 or 2°C; Giudicelli 
et al., 1980; pers. obs.) whereas in other (lower) 
reaches temperature fluctuations will follow the 
air temperature throughout the year (dif¬ 
ferences of some 10 to 15°C between summer 
and winter). 

According to Goodenough (1984), random 
drift and selection pressure allow a “splinter’* 
population to acquire very different allelic fre¬ 
quencies from its parent Mendelian population 
within a very few generations (from about 7 in 
experiments performed by Powell & Richmond 
with Drosophila paulistorum as discussed by 
Goodenough, 1984). Thus, since in the area 
studied environmental conditions vary greatly 
it is very likely that genetic differentiation is due 
to genetic drift, bottlenecking and different 
selective forces acting together. 
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